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There were two objectives in the experiments performed on the 
barotropic model. One was to correct the lag of the troughs and 
ridges noted in the barotropic model by the Naval Weather Research 
Facility, and the other was to introduce a baroclinic development term. 

The approach to these problems was one which could be incorpo- 
rated into the Fleet Numerical Weather Facility model with ease. 

Included herein are the theoretical development of the numerical 
model, a comparison with the FNWF model and charts illustrative of 
the results with various values of the pertinent parameters for 24 
and 48 hours. 

It was observed that the development term causes selective 
development of lows and highs; but when overdone, spurious short 
waves are formed which, however, do not basically change the major 
wave location. 

It was also observed that the function of latitude used here to 
correct the lag in the horizontal advection of vorticity, resulted in 
greatly improved movement of troughs and ridges from 30 degrees 
to 60 degrees of latitude, but could be further refined to produce even 



better results. 
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1. Introduction 



The numerical model for the prognosis of the 500-mb level used 
by the FNWF is based upon the barotropic potential vorticity 
equation Q5] . This model has been carefully '’tuned 11 and engineered 
over a period of several years and is found to produce quite satisfactory 
results. The use of a numerical model derived from a baroclinic 
vorticity equation, however, should improve prognosis by permitting 
development and decay of pressure systems. Due to time limitations, 
it was decided to utilize a quasi-baroclinic model which makes use of 
of the initial ( 1000- 500)-mb thickness and a simplified prognosis 
thereof to compute a Sutcliffe development term in the 500-mb prog- 
nostic equation. No attempt is made to verify the thickness forecasts 
nor the implied 1000-mb prognosis. 

The model developed and tested here is still essentially a one- 
parameter model and constitutes only an attemp to provide for some 
development in the conventional barotropic model, which is the primary 
shortcoming of an otherwise very successful prediction scheme for 
mid- troposphere. 

In addition to the use of a baroclinic term, some experimentation 
was performed on the horizontal vorticity advection term of the 
equation by varying it as a function of latitude. Various degrees of 
smoothing were also applied to the prognosis and the effects on the 
appearance and verification are noted. 
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2 . 



The FNWF model jj] 

The barotropic potential vorticity equation is taken as 



_L Cf*r) 

* 



- o 



( 2 . 1 ) 



where h is the depth of fluid being considered. For application to the 
atmosphere, this layer is considered to be approximately the tropo- 
sphere and to be centered at the 500-mb level. Expanding equation Z. 1 
and multiplying through by h gives 



v-vr-+( hfJ\V-Vh - o , 2 . 2 ) 

The last term is assumed to be composed of three components; 
namely, 

a) a mountain term, 

b) a latitude term, 

c) a temperature term. 

Equation Z. Z then becomes 

-■£((+?)& * V-Vf+V' Vh, 

(i) (a) ID Ott V) (2.3) 

») (?) 

Term (Z) is the Helmholtz term. It is approximated by 




where Z is the height of the 500-mb surface, f is a mean value for (f ), 
h is a mean value for h, and is a constant determined by experiment- 
ation. 



1Z 



Term (5) is a mountain term. The approximation 



v> vtv* -v- v ( 2. si 

n h 

is made where M is the height of the terrain, and M a M is an emperically 
determined constant. With other considerations and simplifications 
this results in 

( V 1 - a) 2|. ^ o.-f M - b£H') 

* ?t "ft ' (2.6) 

where n b M is a tuning constant and H is a (500- 1000)-mb thickness. 

The thickness is carried forward during the forecast a§ an 
independent field with half of the speed of the 500-mb wind. 

In addition to the above, other engineering and climatological 
modifications are made in the FNWF model but will not be discussed 
here. 
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3. The baroclinic model 



The baroclinic model used here is developed from the horizontal 
advective model £ 2 , p. 400] in which temperature changes are assumed 
to result from horizontal advection only. In this model, it is also 
assumed that the vertical variation of the magnitude of the thermal 
wind is the same everywhere and the wind field may be expressed 




(3. 1) 



where ^7 is the vertically integrated mean wind with respect to 



pressure and is the vector difference between the wind at two 



T 



selected isobaric levels. With the geostrophic assumption, is 



the thermal wind. 



By the same assumption 




(3.2) 



Integrating the vorticity equation 




(3.3) 



after substituting from (3. 1) and (3. 2) gives 



i.J, jk(*+wsf)fc+’fcX P)y r +f)ip 

3 £ ^ 1*00 



(3.4) 



The left side leads to 




(3. 5) 
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By the geostrophic approximation 



y = |- kx vl 

f* Lv 

T 

substituting (3.6) into (3. 5) 
leads to 




3.8 becomes 



V l ||-+-J(2,5’)+A‘J(H>f7"H) (3.9) 

The right side of (3.4) leads to a Helmholtz term for which the FNWF 
version is utilized. Introducing tuning coefficients K and K yields 

M J 

to the final operational equation 

<3. ic 



(3.6A) 
(3. 6B) 
(3. 6C) 
(3. 6D) 
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where Z is the 500-mb height, H is the ( 500 - 1000) -mb thickness, and 



A is the FNWF Helmholtz term. Assuming the thickness changes are 
purely advective gives 



If- •* fr'-ftf (3.U) 

where K . is 0. 5. 

4 
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4. The numerical form of the baroclinic equation 



Equations (3, 10) and (3. 11) are solved over a field composed 
of 63x63 grid points which are a distance d apart. The appropriate 
map factor m for the FNWF 63x63 grid is that for a polar stereo- 
graphic projection £l) . In order to utilize the FNWF subroutine 
library and reduce computation time, fixed point arithmetic was used. 
The constants and corresponding scaling factors are: 

d = 381KM = 3. 81 x 10^ cm 
m = 2m 

f = 1. 03125 x 10" 4 sec’ 1 

-2 

g = 980 cm sec 
D = D x 2 1? (D = Z-Z'std) 

H = H x 2 1? 

? =£ x 2’ 9 
x 2 

At = 1 hour = 3600secs 
which result in 



K** t ( H, ^ ) (4. l) 

which is the form of equation (3. 10) used in the FNWF Helmholtz 
subroutine . 

The thickness is carried forward during the forecast as an 
independent field with half of the speed of the 500-mb wind 



§T-*-±V‘VH*‘±T(2Al) 
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(4. 2) 



When expressed in finite difference form, (4. Z) becomes 

& ft, H ) (4.3) 

Both (4. 1) and (4. 3) are computed on an hourly basis. The first 
time step is produced by the method of forward time differences 

D m+1 = D m + AD m (4.4) 

If continuedjthis method would lead to computational instability. For 
this reason, all other time steps are produced by the method of 
centered time differencing 4 

_m+l ^m- 1 ^ k ^_m 

D =D +2AD (4. 5) 

In the event a smoothing operation is performed on a field, a forward 
time difference is used as the subsequent time step. 
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5. Research procedures 

Comparing equations (2. 5) and (3. 10) shows that the FNWF model 
without certain refinements closely approximates the baroclinic model 
when is zero. For this reason, this baroclinic model without the 
development term can be referred to as the quasi-FNWF model; 
however, the FNWF model has additional refinements not included 
here, hence comparisons of verification scores are not strictly valid. 

The data used for the experiments were recent 500-mb and (500- 
1000)-mb thickness fields obtained from FNWF. The actual computing 
was done on the Postgraduate School Computer Facility CDC Model 
1604. 



The results of each experimental run were tested by the pillow 
and RMSE £lj , then analyzed and visually compared with the analysis 



at verification time as well as with other experimental runs. 

X 



Pillow 



RMSE 



4 Cfa**i)* P\ll6*v0 d. 

**' * 



It is apparent that the RMSE, although a rapid method of 
quantitative comparison of fields, will not show a variation of develop- 
ment of systems. It was therefore necessary to use the visual method. 

The experimental procedure consisted of varying the values of K 
or or both. The values of varied from 0 to 2. 0 and also 
included variations of the tuning constants with respect to latitude. 
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In addition to the above, the fields were compared after different 



degrees of smoothing were applied. 
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6. Experimental results 

Consider first the advection term only. The speed of movement 
of the troughs and ridges were modified by the use of as a function 
of latitude. The need for this modification was made evident by a 
study of the FNWF model conducted by the Naval Weather Research 
Facility, Norfolk, Virginia, which indicated consistently slow move- 
ment of troughs and ridges through certain latitudes, a deficiency that 
appears to characterize barotropic models in general and even some 
baroclinic models. 

After several trials, an emperical correction was made for this 
lag by the employment of a function of latitude (SF) as the factor 

.r + 

applied to the horizontal vorticity advection, 

The graph of this function. Fig. 1, shows a maximum of 1. 37 
at 21.4 degrees and a value of one at the equator and pole. 

The effect of the use of this sine function can be seen by comparing 
Fig. 2, the 24-hour prognostic chart without the sine function and 
Fig. 3, in which the sine function is applied. 

There is an obvious displacement of all systems while the 
contour shape remains essentially unchanged. 

To test the effectiveness of this function on independent data, a 
test was run on seven consecutive 24-hour periods from 0000Z 24 March 
to 0000Z 27 March. The displacement of troughs and ridges was 
measured in degrees of latitude and an average deviation from the 
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J?ig. 2 
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Fig. 3 
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TABLE I 



24-Hr Prog 

AVERAGE DISPLACEMENT FROM 
ANALYSIS 

IN DEGREES OF LATITUDE 



LAT 


WITH SINE FUNCTION 
TROUGHS RIDGES COMBINED 


WITHOUT 

TROUGHS 


SINE FUNCTION 
RIDGES COMBINED 


10 


-2. 72(11) 


-2.0(2) 


-2. 62(13) 


-2. 64(11) 


-2.0(2) -2.54(13) 


20 


1. 43(23) 


36(14) 


. 76(37) 


. 61(23) 


-.93(14) .03(37) 


30 


. 70(36) 


-1. 91(23) 


-.32(59) 


-.97(36) 


-2.26(23) -1.47(59) 


40 


. 51(35) 


-2. 40(30) 


-. 83(65) 


-.66(35) 


-3. 57(30t) -2.00(65) 


50 


35(31) 


-2. 15(26) 


-1. 16(57) 


-1. 52(31) 


-3.46(26) -2.40(57) 


60 


-1. 66(18) 


-.95(20) 


-1. 30(38) 


-2. 38(18) 


-2. 15(20) -2. 26(38) 


(# of samples) 
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TABLE II 



24-Hr Prog 

TEST OF SINE FACTOR 
PILLOW* RMSE* 



V. T. 


NO 


SF 


WITH SF 


FNWF 


00 Z 24 Mar 


-1 


93 


-2 


94 


1. 1 


104 


12Z 24 Mar 


-5 


101 


-5 


102 


6. 3 


115. 3 


00Z 25 Mar 


-5 


90 


-5 


91 


6.6 


108. 3 


12Z 25 Mar 


5 


97 


5 


97 


-3. 8 


105. 9 


00 Z 26 Mar 


0 


113 


0 


113 


- .3 


111.1 


12Z 26 Mar 


-2 


117 


-2 


116 


-1. 2 


115. 0 


00Z 27 Mar 


18 


94 


18 


95 


-25. 5 


101. 5 



*feet 
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corresponding analysis was taken for the entire period. 

The results of this test are summarized in Table I. The first 
item observed by this summary was that the movement of troughs 
and ridges was truly in need of an increase at almost all latitudes 
tabulated. It is seen that the use of the sine function produces a 
distinct improvement from 30 to 60 degrees. This is particularly 
evident in the n combined M column which includes troughs and ridges. 

Table I also shows a possibility of too much correction at 20 
degrees . 

The results at 10 degrees are considered to be inconclusive due 
to the small number of samples and the difficulty of defining troughs 
at this latitude. 

Table II compares the pillow and RMSE for the verification dates 
of the test without the sine function, with the sine function and for 
the FNWF product. The use of the sine function produces little change 
in this statistic. 

Turning now to the development of systems, the problem was to 
define a proper term or to determine how much development term 
to use. 

The experimental data was selected at random from the FNWF 
files of height fields and corresponding thickness fields. 

The 24-hour prognostic tests were run from 0000Z 25 Jan 1967, 
Fig. 4, to 0000Z 26 Jan 1967, Fig 5. This date was chosen at random 
and, although not the best of all possible situations to illustrate 
development, does show the effects of the several variations. 
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Fig. 4 
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Fig. 5 
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For purposes of comparison, the following factors are considered: 

a) the pillow(P) and RMSE(R) in feet 

b) illustrative charts 

c) the mean of the deviation of selected low and high central 
values in meters from the analysis, together with the standard 
deviation from this mean. 

It will be seen that a low RMSE does not necessarily indicate a 
good verification and therefore must be carefully used in conjunction 
with other evidence before drawing conclusions. 

For information, the FNWF pillow and RMSE for this date 
were 3. 4 and 115. 9, respectively. The mean of the deviation of lows 
was 7. 8 with a range of 92 to -46 and a standard deviation of 47. 1, 
and the mean of the deviation of the highs was 11.7 with a range of 
31 to -2 and a standard deviation of 14. 0. 

Caution should be used in comparing the FNWF barotropic results 
to the barotropic model used here, where is 0, as this model only 
approximates the FNWF model. However, comparison of cases with 
and without the development term should indicate what effect the 
development term would have on the FNWF model. 

The 500-mb analysis during the forecast period showed little 
significant movement or development with the exception of the low 
pressure system on the northwest coast of the U. S. which moved 
inland and filled, while the low in the Gulf of Alaska deepened by 90 
meters and moved east leaving a trough trailing behind. 
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Case I 



Fig. 6 is the result of using no development term, = 0. The 
general position and intensity of pressure systems are good and 
follow the history intensities. The low on the west coast of the U. S. 
fails to fill while the low in the Gulf of Alaska deepens by only 20 
me ter s. 

mean deviation and range standard deviation from mean 



p 


-7 


lows 


-5. 0(51, -38) 


29.6 


R 


105 


highs 


19. 2(28, 9) 


6. 8 




K 3 


was next 


set at 1. 0. The 


prognosis became more unstable 



with each time step and after 13 hours, the FNWF Helmholtz sub 
routine CQ would not converge. 

Various other values of were tried. It was found that when 
is 0. 5, the model runs for 18 hours, at 0. 3 for 23 hours and 
equal to or less than 0. 2, over 24 hours. 

Since the hourly change of the thickness field is based solely on 
horizontal advection, it is believed that this field becomes less 
representative of the actual field fairly soon and is probably a major 
source of observed instability. For this reason, was divided by 
2 every six hours in the following examples unless otherwise 
indicated. Other periodic reductions of were tried, however, this 
one gave the best results. 

In addition, the FNWF subroutine SAU is employed on the D 

and H fields every three hours to suppress the generation of short 
waves. 
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Pig. 7 
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Fig. 8 
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Case II 



In Fig. 7, is 0. 5. Although the major features remain 
unchanged, compared to Fig. 6, there is some increase in short wave 
development. The mid-latitude lows are from 20 to 40 m deeper and 
highs 10 to 30 m higher. 

mean deviation and range standard deviation 
P -8 lows -8.3(45,-39) 22.7 

R 115 highs 26.4(50, 8) 14.9 

Case III 

In Fig. 8, is increased to 1.0. Further short-wave development 
is noticed in addition to stronger development of both highs and lows. 

The development is selective, as can be seen by the low north of 
Siberia remaining unchanged while the low in Idaho deepens by 50 m 
and the one in the Gulf of Alaska by 80 m. 

mean deviation and range standard deviation 
P -8 lows -16. 8 (2, -43) 15. 7 

R 128 highs 34. 6(62, 8) 22.3 

In the above examples, it becomes evident that the use of the 
development term produces an unacceptable amount of short-wave 
development. FNWF has a subroutine, SAG jflj , which can be used 
to remove some or all perturbations. It is of the form 

A = A + k 7 2 A 

The problem remains as to what size to make k and how many 
times to use the smoother. 
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Pig. 10 
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The smoother was tried with k = 0. 125 and used three successive 



times on the D and H fields every six hours. The results of this heavy 
smoothing in addition to no sine function and at 0 is seen in Fig. 9. 

Case IV 

Many of the systems have been smoothed out of existence and 
others have become unrepresentative of the verifying analysis. This 
example clearly shows that a low RMSE does not necessarily mean a 
good prognosis. 







mean deviation and range 


standard deviation 


p 


-7 


lows 116.0(174, 92) 


32. 1 


R 


105 


highs -29.6(5,-77) 


34. 2 




Case 


V 





Maintaining the heavy smoothing, and setting at 1.0, starts to 
bring the prognosis pattern. Fig. 10, back to the analysis, but not 
enough. 

mean deviation and range standard deviation 
P - 8 lows 110. 3(142, 84) 22.4 

R 107 highs -26. 2(3, -60) 32.0 

Case VI 

In Fig. 11, is 1.0 but is not diminished during the 24 hour 
period. As noted earlier, this would be too unstable to be computed, 
were it not for the use of a heavy smoother. 

Compared to Fig. 9, both highs and lows show much greater 
development, but this is attained with an excessive gain in short 
waves . 
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Pig. 11 
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Fig. 12 
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Pig. 13 
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mean deviation and range 



standard deviation 



P -8 


lows 70. 9(117, 14) 


41. 2 


R 147 


highs 8.4(63, -30) 


30. 5 



Case VII 

Continuing the search for a compromise between development and 
smoothness, is again divided by 2 every 6 hours and the initial 
value of is set at 2.0, which results in Fig. 12. While the central 
heights are closer to the analysis than Fig. 9 or 10, there is still too 
much short-wave development, although less than Fig. 11. 



mean deviation and range standard deviation 



P 


-8 


lows 85. 5 (117, 55) 


23. 9 


R 


128 


highs 6. 2(53, -38) 


33. 0 




Case VIII 






In Fig. 


13, the smoothing coefficient 


is reduced to 1/16 vice 1/8 


while is 


1.0. This is seen to produce development with acceptable 


smoothness 


• 








mean deviation and range 


standard deviation 


p 


-8 


lows 48. 9(69, 28) 


14. 7 


R 


116 


highs 4. 0 (36, - 21) 


18. 8 



Case IX 

Continuing in the same vein, the smoothing coefficient is reduced 
to 1/32 while remained 1.0, Fig. 14. Further development is noted 
but so is considerable short-wave formation. The compromise point 
has been passed. 
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Fig. 14 
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Pig. 15 
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Pig. 16 
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Pig. 17 
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mean deviation and range 



standard deviation 



P -8 


lows . 3(26, - 22) 


17. 5 


R 140 


highs 27. 6(65, 7) 


21.6 



The above findings are now tested for a 48 hour period using a 
verification analysis of 1200Z 25 Jan. 1967, Fig. 15. The FNWF 
pillow and RMSE for this date were -16.0 and 173. 2, respectively. 

The mean deviation, range and standard deviation of the lows were 
-81. 0(119, -230), 119 and for the highs -24. 3(- 10, -39) and 11. 8. 

Case X 

Starting with the heavy smoother: k of 1/8, 3 scans every 6 hours, 
and at 0, results in Fig. 16. This is an extremely smooth chart 
with most baroclinic features missing. A misleading low RMSE is 
again encountered. 

mean deviation and range standard deviation 
P 8 lows 117. 0(165, 81) 102.7 

R 143 highs -154(1 sample) 

Case XI 

Continuing with the heavy smoother and changing to 1.0 produces 
some slight development, Fig. 17, but still far short of the analysis. 

mean deviation and range standard deviation 
P 7 lows 101. 3(147, 66) 108.4 

R 142 highs -154(1 sample) 

Case XII 

Reducing the smoothing factor to 1/16 with of 0 improves the 
chart, Fig. 18, over the previous two results but still not enough. 
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Pig. 18 
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T 




Pig. 19 
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Fig. 20 
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mean deviation and range 



standard deviations 



p 


7 


lows 67. 0( 1 1 2, 37) 


110. 8 


R 


148 


highs -126(1 sample) 


- 




Case 


XIII 






With the smoothing coefficient at 1 / 1 6 and at 1.0, further 


development is observed in Fig. 19, but 


even more would be desi 






mean deviation and range 


standard deviation 


p 


7 


lows 35. 4(84, - 5) 


119 


R 


161 


highs -67. 0(- 24, -110) 


43 



An attempt was made at reducing k to 1/32 while leaving at 1.0. 
This variation produced a 30 hour prognosis before overflow occurred. 
Case XIV 

The final variation, Fig. 20, was the most satisfactory although 
not spectacular. In this variation, k is 1/16 but used twice instead of 
three times every six hours while is 1.0. 







mean deviation and range 


standard deviation 


P 


8 


lows - 55. 4(- 14, - 123) 


136. 5 


R 


168 


highs 5. 0(62, -62) 


52. 4 




The 


results of the above cases are 


summarized in Table III 
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RESULTS OF VARIED DEVELOPMENT TERM VICE SMOOTHING 
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TABLE III (Continued) 
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7. Independent data results 



The last combination of factors, that is, of 1. 0 reduced by one 
half every six hours and a smoothing coefficient of 1/16 applied to the 
height and thickness fields every six hours, was applied to independent 
data of seven consecutive 24-hour prognostic charts. 

These charts covered the period 0000Z 24 March to 0000Z 27 March. 
A comparison was made with full development term of 1. 0, no develop- 
ment term and the FNWF product for these dates. 

Upon summarizing the results of the above tests, it was seen that 
a of 1.0 was probably too high and another group of seven charts 
was produced with a of 0. 7. 

The results of all of the above are listed in Table IV. The pillow 
and RMSE for the variations on the applicable dates is presented in 
Table V. 

Table IV indicates a sharp reduction in the mean deviation of the 
lows from the analysis, however, the choice of a of 0. 7 was too 
much of a reduction from 1 . 0 in the case of the lows but was quite 
good for the highs. A of about 0. 8 would probably have produced 
better results for this period. 

Compared to FNWF, a of 0. 7 results in a good reduction in the 
range of deviations for both highs and lows and a marked improvement 
for the standard deviation from the mean deviation for highs. 
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INDEPENDENT DATA TEST OF DEVELOPMENT TERM 

LOWS (31 cases) HIGHS (12 cases) 

Mean deviation Standard Mean deviation Standard 

from analysis Range deviation from analysis Range deviation 
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TABLE V 



INDEPENDENT DATA TEST OF DEVELOPMENT TERM 
PILLOW AND RMSE (ft. ) 
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8. Conclusions 



The modification of the horizontal advection term by a function of 
latitude as used herein, generally improved the movement of troughs 
and ridges by reducing the lag between prognosis and analysis. 

It appears that a refinement of this function with a displacement 
of the maximum from about ZO degrees to about 45 degrees would be 
more rewarding. 

While compiling data for this study, there appeared to be a 
relation between the zonal index and the lag of the troughs and ridges. 

It is recommended that further study of this phenomena be conducted 
to arrive at a more exact relationship and how it could be applied to 
proper movement of troughs and ridges. 

It is readily seen that, as the name implies, the development term 
produces development in the form of deeper lows, higher highs and 
also, unfortunately, unwanted short waves. 

It results in distinct and selective changes in intensities of systems 
without materially changing the location of these systems. 

It showed marked improvement over the FNWF highs and a lower- 
ing of the range of errors in the lows. 

The inclusion of this term to some degree in an operational 
barotropic model would be beneficial. 

It is suggested that the most desirable approach would be to 
incorporate these empirical modifications in the FNWF 500-mb 
prediction equation while simultaneously carrying out the FNWF 1000-mb 
prediction hour by hour. The resulting H field could then be incorpor- 
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ated into the 500-mb forecasts each hour, which should be far more 



reliable than that produced by simple advection of thickness. 

A development term so produced could be made smaller initially, 
say = 0. 25, but would not be reduced periodically since the thick- 
ness field would presumably stay in phase better with the 500-mb 
height field. * 

This scheme should produce improved long range 500-mb and 
surface prognostic charts. 
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